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Abstract The torsional elasticity of semiflexible polymers like DNA is of 
biological significance. A mathematical treatment of this problem was begun 
by Fuller using the relation between link, twist and writhe, but progress 
has been hindered by the non-local nature of the writhe. This stands in the 
way of an analytic statistical mechanical treatment, which takes into account 
thermal fluctuations, in computing the partition function. In this paper we 
use the well known analogy with the dynamics of tops to show that when 
subjected to stretch and twist, the polymer configurations which dominate 
the partition function admit a local writhe formulation in the spirit of Fuller 
and thus provide an underlying justification for the use of Fuller's "local 
writhe expression" which leads to considerable mathematical simplification 
in solving theoretical models of DNA and elucidating their predictions. Our 
result facilitates comparison of the theoretical models with single molecule 
micromanipulation experiments and computer simulations. 
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1 Introduction 

DNA, which carries the genetic code of living organisms is a semiflexible 
polymer. The elastic properties of DNA are relevant to a number of biolog- 
ical processes [5]. Twist rigidity plays an important role in packaging metres 

J. Samuel 

Raman Research Institute 
Tel.: +91-80-2361 0122 
Fax: +91-80-2361 0492 
E-mail: sam@rri.res.in 



S. Sinha 

Raman Research Institute 



2 



of DNA efficiently in the tiny volume of the cell nucleus, just a few microns 
across. This involves DNA-histone association which makes use of supercoil- 
ing in an essential way. The process of DNA transcription can generate and 
be regulated by supercoiling|31j. 

To understand such effects, there have been single molecule experiments^^ 
which pull and twist DNA molecules to probe their elastic properties. How- 
ever, the interpretation of these experiments has been hindered by a lack of 
understanding of the geometry of writhe, the twisting of a DNA backbone. 
The writhe of a DNA polymer configuration, viewed as a space curve, is a 
nonlocal quantity which has subtle geometric and topological properties. The 
nonlocality of the writhe makes it cumbersome to use in analytical or compu- 
tational work [TTI|331I22[|211I32"I ITIH1. This paper is devoted to understanding 
the writhe for configurations which dominate the partition function: those 
close to the minima of the energy functional. More precisely, we show that 
for the stable solutions of the Euler-Lagrange equations (local minima of the 
energy) the writhe can be expressed as a local integral. This result can be 
used to develop an analytical approach to the computation of the partition 
function. 

The most popular theoretical model for describing semiflexible polymers 
is the worm-like-chain (WLC) [TO] . The WLC is the " harmonic oscillator" 
in the field of semiflexible polymers: it is a simple model which deftly cap- 
tures much of the physics. Given the importance of semiflexible polymers 
to biology, from lipids to cytoskeletons to DNA, the WLC deserves to be 
much better understood than it presently is. Two extreme situations are 
relatively well understood. If the polymer is highly stretched and nearly 
straight, one can use perturbation theory about the straight line to calcu- 
late its elastic properties [1 9 ,28,29, 2"SlfT5] . Treatments also exist [33] in the 
opposite extreme when the polymer is wrung so hard that it buckles and 
forms plectonemic structures [TB], which are stabilised by the finite thick- 
ness of the polymer. The transitional regime where the polymer is neither 
straight [28 , 19] nor plectonemic [T5] is the subject of this paper. Perturba- 
tion theory about the straight line is not applicable and the simplifications 
that arise from the energy dominated plectonemic regime do not obtain. We 
address this regime by finding the minima of the energy functional in the 
partition function. These minima come in two families- the straight line fam- 
ily and the "writhing" family, which deviates considerably from the straight 
line. We study these configurations which dominate the partition function 
in the transitional regime. We show that these curves admit a local writhe 
formulation. This leads to a considerable simplification of the theory. 

The paper is organized as follows. Sec II presents a recapitulation of 
the geometrical notion of writhe of a simple space curve and how it can be 
applied to understand DNA elasticity. Sec III deals with the mechanics of 
semiflexible polymers using the classical mechanics of a top [13]. We prove our 
main result here: the polymer configurations which dominate the partition 
function admit a local writhe formulation. Brief remarks regarding perturba- 
tions about the saddle point and stability are made in Sec IV. We end with 
some concluding remarks in Sec. V. 
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2 Geometry 

In this section we deal with the essential geometric notions that determine 
the twist elastic properties of linear molecules like DNA. This section is 
a recapitulation of some material that has been presented earlier |25j and 
is being summarised here for the reader's convenience. The experiments of 
Stride et al [31] are done on an open segment of DNA which is attached to a 
glass slide at one end and a magnetic bead on the other. The bead is pulled to 
stretch the molecule and turned to "wind it up" and induce supercoiling. For 
mathematical convenience, we close the open segment of DNA by a reference 
ribbon (a ribbon is a framed space curve) in a manner described in |25j (see 
Figs. 4 and 5 of Ref[Tl] and references therein). The fixed reference loop 
(called the passive part) is held fixed and our discussion concerns only the 
active part of the curve. (It is also possible to choose the reference ribbon as 
coming in from infinity along the z direction and going off to infinity along 
the z direction. The precise choice of reference does not matter.) 

The basis of the analysis of DNA supercoiling is the celebrated relation 
[SEEHI]: 

Lk = Tw + Wr (1) 

which relates the applied link (the number of times the bead is turned in the 
experiment) to the twist (the twisting of the polymer about its tangent vec- 
tor) and the writhe (the bending of the tangent vector) . The problem neatly 
splits [TT1I25] into two parts: a relatively straightforward (local) description 
of the twist elasticity, and the somewhat harder (nonlocal) treatment of the 
writhe. 

The writhe is a non-local quantity defined on closed simple curves: Let 
the arc length parameter s range over the entire length Lq of the closed 
ribbon (real ribbon + reference ribbon) and let us consider the curve x(s) to 
be a periodic function of s with period Lq. Let R(s,er) = x(s + a) — x(s). 
We write i{s) = ^ for the unit tangent vector to the curve. Since the curve 
x(s) is simple, R(s, a) is non- vanishing for a ^ 0, Lq and the unit vector 
R(s, a) is well-defined. It is easily checked that as a — > {0, L }, R — > {i, —i} 
respectively. The Calugareanu- White writhe is given bv [7, 6, 34 TO, XT) 



1 I , f , r dR(s,o-) dR(s,a), - 
W C w = — & ds da[ — y—^ x \ > ] ■ R. (2) 
4tt J J 0+ ds da 

The non-locality of Wcw makes it difficult to handle analytically. However, 
the key point to note is that variations in Wcw are local [111110) . Let C(a) be 
a family of polymer configurations parametrised by a of simple closed curves 
with writhe Wcw{ a )- Taking the a derivative of Eq. ^ (which now depends 
parametrically on a) we find that the resulting terms can be rearranged to 
give 

dWcw 1 f Lo , ,dt dt t 



da 2tt Jy ds da 



which clearly has the interpretation of the rate at which t(s, a) sweeps out a 
solid angle in the space of directions. Note that Eq.Q which gives the change 
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in writhe is a single integral ;3 , 8 and therefore a local, additive quantity. This 
can be expressed as a sum of the changes in the active parts of the curve and 
therefore, since the passive part is unchanged in the variation, the change in 
writhe is entirely due to the active part < s < L, where L is the contour 
length of the active part of the polymer. Eq. (|3]) can be rewritten as: 



where Q is the solid angle enclosed by the oriented curve {i(s) |0 < s < L} on 
the unit sphere of tangent directions as s goes from to L [18 . Note that Q 
is only defined modulo An: for a solid angle ft to the left of the oriented curve 
t(s) is equivalent to a solid angle (An — f2) to the right, dfl/da is, however, 
well defined and local. Integrating Eq. (Q| we arrive at [TT1ITU] : 



where n is an arbitrary integer. The constant of integration is fixed by noting 
that the writhe of a simple planar curve vanishes (EqJ2]). 

There is a quantity one can construct from the writhe which is well defined 
on all curves (not just simple ones) 



is a complex number of modulus unity which we call the wreathe. When a 
curve is passed through itself, Wcw jumps by 2, but the wreathe is un- 
changed. We can therefore smoothly extend w to all closed curves including 
non-simple curves. Unlike the writhe, the wreathe w is a local quantity. The 
wreathe can be exploited to define a new quantity, the Fuller "writhe" 
HOj for a class of curves. (Readers familiar with the geometric phase in quan- 
tum mechanics [2"l l2"Tll2"51l2~4"] will recognise the relation between the wreathe, 
the Fuller writhe, the Berry phase and the Pancharatnam connection.) Let 
us fix a direction (say —z, which we call "south") and define "south avoid- 
ing" curves to be those for which the tangent vector never points south Q. 
The Fuller "writhe" is defined for south avoiding curves. Let us define the 
"wreathe angular velocity" A a = —iw^ 1 ^- = l/2dQ/da on the family of 
curves C a parametrised by a. Let us choose a fiducial curve C*, of length L, 
whose tangent vector is identically north pointing. We take this as a standard 
for the "active part" of the curve, the "passive part" , of course being held 
fixed in the discussion. Observe that all south avoiding curves are deformable 
to the fiducial curve C*. One simply deforms the tangent vector t(s) along 
the unique shorter geodesic connecting t(s) to the north pole. We now define 
the "Fuller writhe" as 




(4) 



Wcw(a) = —H-l + 2n 



(5) 



w(C) — exp[inWcw (C)} 



exp[iS7/2] 




daA a — L 



1 The sphere of tangent directions has perfect rotational symmetry. The choice 
of the south direction — z is motivated by the fact that in a real experiment one 
applies force and torque in a specific direction which we choose to be along +z. 
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Writing the unit tangent vector as i — (sin 9 cos 4>, sin 9 sin </>, cos 9) , J Q da ^ 

can be written as J ds^j(l — cos 9) for all curves for which the tangent vector 
never points towards the south pole of the sphere of tangent directions. We 
can therefore write [31[S] a local writhe formula on such curves which we call 
"south avoiding curves" : 

W F = ^~ [ds(l- cos 0)^-1. (6) 

Z7T J as 

While (Eq. [5]) is expressed in local co-ordinates on the sphere, it has a clear 
geometric meaning: 2ir(l + Wf) is equal to the solid angle swept out by the 
unique shorter geodesic connecting the tangent vector t to the north pole. 
This definition is explicitly not rotationally invariant, since it uses a fixed 
fiducial curve C* and singles out a preferred direction. 

To summarise, Wcw is defined on all self-avoiding curves, Wf on all 
south avoiding curves. When a curve passes through itself, Wcw jumps by 
two and when the tangent vector to a curve swings through the south pole, 
Wf jumps by two units. The writhe is a real number which has both geometric 
and topological information. The topological part is the integer part of w <? w 
and the geometric part is the fractional part of w< ^ w . The geometric part 
is completely captured by wreathe but the topology is lost since wreathe is 
insensitive to changes in writhe by 2 units. From the definitions it is clear 
that on a family of curves C a parametrised by a, 

w" 1 dw dWcw dW F 

— i = = (7) 

7r da da da 

So, continuous changes in writhe are the same whether measured by Wcw, Wf 
or w. We have adjusted the reference ribbon and the definition of Fuller 
writhe so that these two notions of writhe Wcw and Wf agree on the ref- 
erence curve. Note that w is well-defined on all curves. For a closed circuit 
of curves {C a , (0 < a < l),Co = Ci}, we have the "closed circuit theorem" 
|25j : In a closed circuit of curves, the total number of south crossings and 
self crossings (both counted with sign) are the same. This follows easily by 
integrating (Eq. [7]) along a closed circuit. 

This definition of "Fuller writhe" is motivated by a Theorem of Fuller [TT1 
110] . which uses a reference curve and states precise conditions (which are 
frequently misunderstood [20 ) under which the writhe difference between the 
two curves can be written as a local integral. Under deformations of the 
reference curve which are both south avoiding and self avoiding (we follow 
[2T3] in calling these "good deformations"), Eq. (0 implies that the equality 
Wcw = Wf is maintained. These deformations are the ones which satisfy 
the conditions of Fuller's theorem [Eq. 6.4 in Ref. [H]]. We refer to curves 
which can be reached from the reference curve by "good deformations" as 
"good curves" . 

To avoid misunderstanding, we remark that the "Fuller writhe" is a dis- 
tinct quantity from Wcw and is not equal to the Calugareanu- White writhe. 
This is obvious since the first is local, while the second is not. Strictly speak- 
ing Fuller only defines a local writhe on "good curves" , where it is equal 
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to the CW writhe. We find it convenient [25] to extend Fuller's defintion to 
all south avoiding curves and not just "good curves". On "good curves", 
the two notions of writhe agree and we have Wf equal to WcWi but the 
two quantities are different in general. The main claim of [25] is that "good 
curves" dominate the partition function in a region of the parameter space. 
We will substantiate this claim in the next section by identifying the domi- 
nant configurations that contribute to the partition function in a saddle point 
approximation and showing that they are "good curves" . 

In the high force region [2"S11T9] this point is obvious since the polymer 
fluctuates about the straight line configuration which is a good curve. The 
content of the following section is that even at lower forces (or at large 
torques, close to buckling) when the polymer is far from straight, the minima 
of the energy are "good curves". The analysis involves a combination of 
mechanics, geometry and topology. 



3 Mechanics 

Scmiflcxible biopolymers like DNA are subject to thermal fluctuations in a 
cellular environment. Therefore the natural theoretical framework for study- 
ing the elastic properties of such systems is statistical mechanics, which in- 
volves a competition between energy and entropy. In the experiment, one 
controls the force F on the bead and the link Lk, the number of times the 
bead is turned. This brings up the theoretical problem of computing the par- 
tition function Z(F,Lk), the number of configurations (counted with Boltz- 
mann weight) of the ribbon which have a given link (the link distribution). 
We know that the link distribution is a convolution of the writhe distribu- 
tion and the twist distribution 3 , 28 . The latter is a simple Gaussian integral, 
which is easily computed. The real problem is to compute the writhe distri- 
bution of an inextensible curve. By simple transformations |25j the problem 
of determining the link distribution of an inextensible ribbon can be reduced 
to computing the writhe distribution of a space curve. The partition function 
to be computed is 

Z(F, W)=J2 exp[-S(C)/k B T]S(Wcw - W) (8) 
c 

where F is the force, £{C) the energy (defined below in eq.©) and Wcw 
is the writhe. The sum is over all allowed configurations of the polymer[25 . 
We expect the sum to be dominated by configurations which minimise the 
energy functional. Configurations near the minimum will also contribute due 
to thermal fluctuations. We first do a purely classical elastic analysis taking 
only the energy into account, while ignoring the entropy[10, 15,5,8 . This al- 
lows us to draw on physical intuition derived from the elasticity of beams, 
cables, telephone cords and ribbons and paves the way for a fuller treatment 
which incorporates thermal fluctuations around the classical solutions. We 
analyze the classical elasticity of a torsionally constrained stretched semiflex- 
ible polymer. 
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The problem can be reduced to minimizing the energy (we use a dot for 
the s derivative, it reinforces the analogy [14j with the dynamics of tops) 



£[i( s )} = o / l - tds ~ / F - ids ( 9 ) 



2 .10 Jo 



of a space curve x(s) whose tangent vector is t(s) = subject to a writhe 
constraint 

W CW = W (10) 

where, Wow, the Calugareanu- White writhe has been defined earlier. The 
tangent vector t(s) is varied subject to the boundary conditions t(0) = t(L) = 
z fixing the tangent vector to the curve at both ends. As mentioned clearly in 
[TS1I2B1I2S] j the full problem involving the link distribution (with finite twist 
elasticity) can be reduced to computing the writhe distribution, which is 
essentially the limit of infinite twist rigidity (C — ► oo). 
Using the method of Lagrange multipliers we arrive at 

8£(C) = 2itt6Wcw (11) 

where 2itt is a Lagrange multiplier, r can be physically interpreted as a 
torque. The variation in the energy 5£ is given by 

8£ = - [ {i+F)-St (12) 
Jo 

The variation of Wcw ( a s explained earlier (j3])) is a local quantity 

SWcw = f [ixh- Sids (13) 
27r Jo 

and the Euler-Lagrange (E-L) equations are 

- i - F = r(t x t) - 7*, (14) 

where the term jt arises since St ■ t — 0. Note that the E-L equations are 
local differential equations rather than the nonlocal integral equations that 
one may have expected from the variation of a non local quantity Wcw- 
Thus, the variations of writhe W are local. As we will see, this is the key idea 
that allows us to prove that "good curves" dominate the partition function. 

The statics of a twisted beam (or cable or DNA) is formally similar to 
the dynamics of a heavy symmetrical top, a fact that has been well known 
since Kirchoff |14| . The analogy is useful for integrating the Euler-Lagrange 
equations. We use quotes for the analogous top quantities. The "kinetic en- 
ergy" is given by T = \t-t and the "potential energy" is V — F ■ t. The total 
"energy" 

H = T + V = he 2 + sin 2 9<j) 2 ) + F cos 8 (15) 
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is a "constant of the motion" as is the z component of the "angular momen- 
tum" 

J z = [t x i) z - ri z = sin 2 0(f) — t cos 9, (16) 

where we have introduced the usual polar coordinates on the space of tan- 
gent vectors t — (sinScos^jSinflsini?!), cos 8). Using these "constants of the 
motion" we reduce the problem to quadratures as described in |13| . The basic 
equations are 



^ 2 sin 6 



■ j z + Tcose 

<P= — ■ (18) 

sin 9 

Setting u — cos 9, we find that 

u 2 = f(u) (19) 

where f(u) is a cubic polynomial in u. Our boundary conditions imply that 
9 = is a point on the solution. Since 9 has to be finite at 9 = 0, we have 
J z = —t. The form of f(u) simplifies to 

f[u) = (1 - u)[2(H - Fu)(l + u)-t 2 {\- «)]. (20) 

(US]) can be rewritten as 

(21) 



1 + u 



From (fT7|) we also find that T-L > F, else 9 2 would have to be negative at 
9 = 0. We see that f'{u)\ u=1 = -A{U - F) < at it = 1. At it = -1, 
/(— 1) = — 4t 2 < and so f(u) has one physical root at u = 1 and the other 
at Mo between —1 and 1: —1 < uq < 1. The third root ui is generally (except 
for iti = ±1) outside the physical range and U\ > 1 for positive force and 
Ui < —1 for negative force. Since f(uo) — 0, we find from (|2"01) the useful 
relation 

t2 = 2(H-Fu )(l+u ) 

1-UQ [ ' 

The "motion" is confined between the turning points 1 and uq'- uq < u < 1. 
We write ito = cos^o, where is the turning point for 9. Just as in the top 
problem, we find that the solution u(s) is periodic with a period given by 

P = 2 f 4t=. (23) 

Because of the boundary conditions, we have L = nP, where L is the length 
of the polymer and n is an integer. L and n fix the integration constant %. 
We will see in the next section that solutions with n > 1 are not minima of 
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the energy. We therefore restrict ourselves to n = 1. Our interest is solely in 
stable solutions, that is local minima of the energy. 

General solutions to these equations can be found in terms of elliptic func- 
tions but a pedestrian approach gives more immediate insight. The simplest 
solution is the straight line u(s) = 1 for all s, which solves the E-L equations 
(fT4"| for all F, t with F — H for any L. However, the straight line cannot ac- 
commodate writhe and we have to look at the "writhing family" of solutions 
parametrised by F, r. These start from u(0) = 1 at s = 0, and have a turning 
point at s = sq (u(sq) — uq) and return to u = 1 at P = 2sq = L. These 
solutions can be explicitly written down in terms of elliptic functions and 
vary continuously with r in the allowed range |r| < yj AF. Integrating (Eqs. 
I17l2ip gives (u(s), 4>(s)) which gives us the solution to the E-L equations t(s). 

From the "top" point of view, t(s) is the solution. However, in polymer 
physics, t(s) is only a convenient way of describing the polymer configuration 
x(s), which is given by 



The use of t(s) rather than x(s), in the variational problem (Eq. O leads 
to a considerable simplification, since the equations of motion are second 
order rather than fourth order. However, there is a price to be paid. After 
solving for the "tantrix" t(s) (plotted in a special case in Fig. 1), we still 
have to integrate (Eg. [Ml) to produce the actual polymer configuration x(s), 
(Fig. 2) in real space. In the polymer context it is meaningful to ask: does 
this configuration intersect itself? Such a question is not significant for tops. 
Nevertheless, as we will see, the motion of a heavy symmetrical top helps us 
answer this question. 

We prove the main result of this paper: all stable solutions of the E-L 
equations are "good curves" in the sense of the last section. For the straight 
line family, the result is obviously true since the straight line is the reference 
curve. The effort here is showing that all stable members of the writhing 
family t(s) parametrised by (i* 1 , r) are "good curves" . To do this we prove 
that the writhing family is nowhere self intersecting or south pointing for 
— I < uq < I. Since u = 1 is the straight line t(s) — z, which we take as 
the reference curve, this proves that the writhing family can be deformed 
(via a deformation of the parameter uq) to the reference curve by "good" 
deformations and the result follows. Below, we explicitly exclude the case 
uq = — 1, which will be treated separately in the next section. 

It is evident that the writhing family is nowhere south pointing for — I < 
uq < 1 since the turning point 9q < -k and the tangent vector t(s) never 
points to the south pole. To prove that the writhing family is nowhere self 
intersecting requires considerably more work. The proof consists of finding 
a direction V along which the tangent vector t to the polymer has positive 
component. This implies that the polymer never bends back to intersect itself. 
The proof uses the E-L equations in an essential way but does not rely on 
the explicit form of the solution. 

Wc first orient our axes by rotating about the z axis so that 4>(sq) = and 
the tangent vector t(s) = (sinO(s) cos 4>(s), sin 9 (s) sin 4>(s), cos 0(s)) lies in 
the x—z plane at the turning point. Consider for each uq = cos 9q the constant 




(24) 
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Fig. 1 The black curve (tangent indicatrix) traced on the unit sphere (green online) 
by the unit tangent vector to a long polymer for a force F = .11 and torque r = .4. 
The dots (in red online) indicate the tips of the unit vectors V (marked V) and z 
(marked N for North) . 




Fig. 2 Figure shows the polymer configuration in real space for the same param- 
eters as Fig. 1. The curve has been thickened, shaded and truncated for easy three 
dimensional visualisation. 



11 



(s independent) unit vector V = (— cos#o> 0, sin#o)- and the function 

g(s) = t(s) • V. 

Note that the function g(s) is symmetric^ about s — so within the period 
< s < 2sq = P = L. We show that a) either g(s) vanishes identically or 
b) g(s) > for all s ^ s . We need only consider the range s < s < 2s , 
so that 9(s) > 0,u(s) < 1. g(2so) is fixed by boundary conditions to be 
sin 0q > 0. The positivity of g(s) immediately implies our central result (see 
below): The polymer configurations under stretch and torque are such that 
the tangent vectors to the polymer all lie in one hemisphere of the sphere of 
tangent directions. 

The proof proceeds in several steps. 
Step 1: g(s) vanishes at the turning point so and has a local minimum at Sq 
if it does not vanish identically. 

Proof: Remembering that 6 = at the turning point, we find 



and 

from which follows 



t(s ) = (sin d ,0, cos 6> ) (25) 
't(so) = (0, sin 6 4>,0) (26) 

g(s ) = *>o)-V = (27) 
9{s ) = £(s ).V = (28) 

So So is a stationary point of g(s). To show that g(s) has a local minimum, 
we compute g(s) using the E-L equations (|T4|) and (|21I25I26I) to find 

g( S o)=sm9 (-^-F) (29) 
1 + u 

Substituting for r 2 (I22p and simplifying gives 

g( S o) = ^^(m-F(l + u )) >0 (30) 

since W > F. The equality occurs only if sin^o vanishes (uq — ±1) in which 
case g(s) vanishes identically, uq = 1 describes a straight line segment, which 
is the reference curve and therefore a "good curve" . uq = — 1 has been ex- 
cluded from consideration. Below we will exclude both cases (uq = ±1) in 
which g(s) vanishes identically. This proves that sq is a local minimum of 
g(s) for (u ^ ±1). 

Step2: At stationary points of g(s), 

g(s)u(s) > (31) 

for S Sq. 



2 This follows s ince 8(s) is an even function of s — so and 0(s) is an odd function 
of s — so- (see eq. (|21|l 'l. Briefly, the orbit is symmetric about its turning point sq. 
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Proof of Step 2: Since g(s) — i. V = at a stationary point, t is orthogonal 

to both V and t and can be written, t = AV x t for some nonzero A. "Angular 
momentum" and "energy" conservation give us 



i.t = 2{H - Fu) 



(32) 



(t x t) z 



t{1-u), 



(33) 



which can be rewritten as 



l^| 2 (l-ff 2 ) 



2(H- Fu) 



(34) 



and 



A(sin 8q — gu) 



t(1 - u) 



(35) 



In Eq. (|34[) both sides are positive since u < 1. Squaring (|35p and eliminating 
\A\ 2 and r 2 between (|22I34I35[) . we find that at a stationary point of g(s), 
g(s) and u(s) satisfy 



where X(u) = (1 — u) 2 /(H — Fu) is a strictly decreasing function of u in the 
range of interest. From X(u) < X(uq) for u > uq it follows that 



which implies (|31[) . 

StepS: so is a global minimum of g(s): g(s) > for s ^ sq. 

Proof: note that since So is a local minimum of g(s) (from Step 1), as s 
increases from sq, g(s) increases from (since s = sq is a minimum). If there 
are no more stationary points g(s) remains positive and we are through. If 
the next stationary point occurs at si > sq, g(s)\ s=Sl — we have g{si) > 0. 
It follows from Step 2 that u(si) > 0. Since u(s) is a monotonic function of 
s in the range (so < s < 2sq), we have u(s) > for s > s± in this range. 
All minima of g(s) for s > si, therefore have g(s) > (again using Step 2). 
It follows that g(s) > for all s (except the one point sq, where g(so) = 0). 
From the symmetry of g(s) it follows that g(s) is positive everywhere except 
at the turning point (so, where it vanishes). 

This immediately implies our main result since if any of the minimum 
energy configurations intersect themselves, x(si) = x(s2) for S2 > s%, we 
have 



(sin6» -gu) 2 < 



(H - Fu ){l + u ) 
(1 - u ) 



X(u) 



(36) 



(sin 6 Q - gu) 2 < (sinflo) 2 



(37) 




contracting this equation with V gives us 
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which is a contradiction since g(s) > for all s ^ so- We have thus proved 
that the writhing family does not intersect itself for —1 < uq < 1. Since uq 
can be continuously deformed to 1 along the writhing family, we conclude that 
the family consists of "good curves" which admit a local writhe formulation. 
We now consider variations around the solutions of the E-L equations to test 
for stability. 



4 Fluctuations around the saddle point 

In a saddle point approximation to the partition function the dominant ther- 
modynamic contributions are expected to arise from the global minimum 
of the energy. We can also expect contributions from local minima due to 
"metastable" states, when the time scale of the experiment is not large 
enough to find the true minimum. All minima, local and global will satisfy 
the E-L equations and so will be among the solutions described in the last 
section. However, not all solutions to the E-L equations are minima of the 
energy. There are some configurations which, if perturbed, will go to nearby 
states of lower energy and decay. These configurations do not contribute sig- 
nificantly to the partition function. To weed these out we need to consider 
the stability of solutions by perturbing the classical solutions to see if there 
are neighboring lower energy configurations. (We do not concern ourselves 
here with zero modes arising from symmetry breaking. They contribute an 
innocuous multiplicative constant to the partition function.) 

For example, we note that solutions with L = nP and n > 1 are unstable. 
This is seen by explicitly constructing neighboring configurations of lower 
energy. Suppose n > 1 and consider the configuration {u id(s)) f \d(s)} which 
solves the Euler Lagrange equations. Construct a new configuration 



where a may be an arbitrarily small real number. The new configuration 

1. satisfies the boundary conditions, 

2. has the same energy as the old one (since the energy is a local integral), 

3. has the same writhe as the old one (since changes of writhe are local 
integrals), 

4. is not a solution to the £L equations at s — P (because of the break in 
the curve at this point). 

It follows that there are neighboring configurations with lower energy than 
the old configuration and the same writhe. We conclude that the solutions 
with n > 1 are not local minima of the energy. A similar remark applies to 
the straight line family if the torque exceeds a critical torque r c given by 
(t c /2) 2 = F + (f-) 2 , the classical condition for torsional buckling of a rod. 
This can be directly seen by perturbing the straight line. For r > r c there are 




(38) 



Uncw(s) =«old(s) 
</?ncwO) = lfio\d(s) + a 



for s > P 



(39) 
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no stable solutions of the Euler Lagrange equations and the polymer must 
double back on itself and enter the plectonemic regime. 

We now turn to a discussion of the exceptional member of the writhing 
family (uq = —1) characterised by r = 0. Note that as r ranges from — r c 
to r c , the writhing family goes in a closed circuit from the straight line to 
the straight line. The closed circuit has one south crossing (at r = or 
equivalently uq — —1). From the "closed circuit theorem" [23], it follows 
that the circuit must also have one self crossing. Since none of the other 
members of the writhing family have self intersections, it follows that self 
crossing must occur at r = 0. Thus the closed circuit has one self crossing 
and one south crossing, both occuring at the same value of the parameter 
r = 0. The writhing family curve r = does not belong to the configuration 
space of either self avoiding or south avoiding models and both Wcw an d 
Wf are ill defined. One can expand the configuration space to allow all curves 
and use a formulation in which energy is minimised for fixed wreathe. It can 
then be shown by an argument very similar to that given above p8[) that 
the curve is unstable. By rotating the second half of the curve by a small 
angle a relative to the first, keeping the energy and wreathe the same, we 
can produce neighboring curves which are not solutions of the E-L equations 
and can therefore be perturbed to lower their energy. 

The writhing family curves in the regime < |t| < r c , include stable 
solutions, and we can take into account thermal fluctuations in a Gaussian 
approximation about the local minimum of energy. Since any stable curve 
of the writhing family is at a finite distance and energy from south and self 
intersection, we can expect small fluctuations around the original curve to 
be " good curves" . Large fluctuations may lead to "bad curves" , but their 
contributions to the partition function will be exponentially suppressed by 
the Boltzmann weight, as made clear by the example of |26j . and not affect 
the statistical predictions of the model. Thus the writhing family can be 
used as a base for a clean discussion of thermal fluctuations without getting 
embroiled in the topological subtleties of writhe. Such a treatment goes far 
beyond perturbation theory around a straight line since the writhing family 
is far from straight. 

To summarize, our analysis shows that there is a range of force and torque 
relevant to single molecule experiments for which the polymer configuration 
curves are "good" curves. 

5 Conclusion 

In this paper we present a classical mechanical approach to biopolymer elas- 
ticity based on well known analogies between twisted polymers and tops. 
Our main result is that the dominant configurations in the partition func- 
tion are "good curves" in a topological sense. This result permits a local 
formulation of the writhe and a non perturbative attack on the problem of 
twist elasticity by incorporating thermal fluctuations around the writhing 
family. Such an approach is expected to work well in the energy dominated 
regime, when the polymer is twisted and stretched. In this paper, we have 
restricted ourselves to the mathematical and conceptual issues regarding the 
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use of writhe. Our main result is useful in further developing the theory [30]. 
Applying these ideas to the WLC leads to detailed predictions about the 
statistical mechanical behaviour of twisted stretched polymers. These will be 
described elsewhere. 

The experiments of Strick et al [3T] have been theoretically analysed by 
Nelson [T5] and Bouchiat and Mezard [3J . The analysis of Nelson is perturba- 
tive about the high force limit (see also [55]) and not controversial. However, 
Bouchiat and Mezard [3] ventured into the non-perturbative regime of low 
forces. By using a local writhe formula they were able to achieve agreement 
with experiments. However, this step is open to criticising 2 . In fact, the 
"derivation" [3] of a "local writhe formula" is flawed since it is based on 
the wrong assumption that the Euler angles are continuous. In [25] . it was 
suggested that in an approximation, the success of Bouchiat and Mezard's 
model could be due to the fact that the partition function was dominated 
by "good curves". The present work establishes unambiguously that this is 
indeed the case. In [26] we have argued that, in the high Link regime, just 
as a self-avoiding polymer winds around itself forming plectonemes, form- 
ing "australonemes" . The energies of these two configurations are similar, 
although the configurations are not. As a result we find that 

Z F (r, Lk) « Z cw {r, Lk) (40) 

i.e. the partition functions based on local and non-local formulations of writhe 
are approximately equal to each other for a wider range of parameters than 
one may have naively expected. The parameter space consists of (F, r, L) 
and the energy dominated regime includes short polymers (L of the order of 
persistence length), stretched polymers(-F somewhat greater than t 2 /4) and 
wrung (or plectonemic) polymers (t 2 /4 somewhat greater than F). 

We expect this work to be of interest to researchers in this interdisci- 
plinary field of biopolymer elasticity. More work is necessary to investigate 
the transition from the energy dominated regime to the entropy dominated 
regime. Analytic methods, numerical experiments and real experiments, all 
have a role to play in this endeavour. Our mathematical analysis of this prob- 
lem relies on insights gained by Fuller [111110] and uses them to understand 
the entropic elasticity of twisted polymers. 
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